arXiv: 1507.03813v2 [cond-mat.mes-hall] 4 Sep 2015 


Infrared absorption of closely-aligned heterostructures of monolayer and bilayer 

graphene with hexagonal boron nitride 


D. S. L. Abergel^’^ and Marcin Mucha-Kruczynski^ 

^Nordita, KTH Royal Institute of Technology and Stockholm University, 
Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden 
^Center for Quantum Materials, KTH and Nordita, 

Roslagstullsbacken 11, SE-106 91 Stockholm, Sweden 
^Department of Physics, University of Bath, Claverton Down, BA2 3EL, United Kingdom 

We model optical absorption of monolayer and bilayer graphene on hexagonal boron nitride for 
the case of closely-aligned crystal lattices. We show that perturbations with different spatial sym¬ 
metry can lead to similar absorption spectra. We suggest that a study of the absorption spectra as 
a function of the doping for almost completely full first miniband is necessary to extract meaningful 
information about the moire characteristics from optical absorption measurements and to distin¬ 
guish between various theoretical proposals for the physically realistic interaction. Also, for bilayer 
graphene, the ability to compare spectra for the opposite signs of electric-field-induced interlayer 
asymmetry might provide additional information about the moire parameters. 

PACS numbers: 78.67.Wj,73.22.Pr,73.21.Cd 


Graphene and hexagonal boron nitride (/i-BN) are two 
representatives of a new c lass of materials, the two- 
dimensional atomic crystals.^^^ Both are made of atoms 
arranged in a single layer of hexagons. However, while 
graphene, made of carbon, is a gapless semiconductor 
with a linear dispersion relation for electrons in the vicin¬ 
ity of the two inequivalent corners of the hexagonal Bril- 
louin zone,l^ /i-BN is an insulator with a band gap of 
6eV.® It is now possible to place graphene flakes on /i-BN 
substrate^ or build stacks with alternating graphene and 
/i-BN crystals.l^ Because of the small difference in lattice 
constants, heterostructures with closely-aligned crystal 
structures of the two materials, form a moire pattern in 
STM imaged and the dispersion of graphene electrons is 
reconstructed into minibands as a result of the periodic 
moire perturbation.l^HISl 

Despite con siderable theoreticaP^^^ and 
experimentaP2l2ll effort, little is known about the 
actual van der Waals interaction between graphene and 
/i-BN. While considering only the inversion-symmetric 
scalar potential was enough to explain the ^alitative 
features of the initial STM measurements,® a more 
complex model involving three inversion-symmetric 
term^l was used to mod el so me of the transport and 
capacitance measurements.®^ The same approach was 
used to demonstrate the importance of the lattice¬ 
mixing perturbation term to explain optical absorption 
measurements.!^ However, recent observation of topo¬ 
logical currents in graphene//i-BN heterostructur^^ 
means that the moire perturbation contains a significant 
inversion-asymmetric mass ternP^ while recent theoret¬ 
ical papers suggest that the inversion-asymmetric part 
of the moire perturbatio n migh t be comparable to the 
inversion-s ymmetric p artl ^^ l ^^ l ^^^ ^H article, we use 

four modelJi®i^E^!^ developed to describe the effect of 
the /i-BN substrate on graphene electrons to investigate 
the infrared optical absorption of monolayer (MLG) and 
bilayer (BLG) graphene on /i-BN. We show that phys¬ 


ically different moire perturbations can lead to similar 
absorption spectra, limiting the amount of information 
about the moire perturbation that can be extracted 
from the spectrum. However, we suggest that careful 
examination of the absorption spectra for doping in the 
vicinity of four holes per moire unit cell will provide 
insight into the character of the moire perturbation and 
an experimental technique for distinguishing between the 
various theoretical proposals for the moire perturbation. 
Moreover, we show that for BLG on /i-BN (where, in 
contrast to MLG on /i-BN, the optic al absorption has 
not been investigated previousljfflEIl) ^ it is important 
to consider the trigonal warping of the electronic band 
structure when interpreting the absorption spectra, as it 
leads to shifts of the miniband edges and hence of some 
of the spectral features. We also discuss the BLG//i-BN 
absorption spectra in the presence of an interlayer 
asymmetry due to a perpendicular electric field. 

For a MLG/BLG crystal with a lattice constant qq = 
2.46A placed on top of a /i-BN flake with a lattice con¬ 
stant a/i_BN = (l + i5)aG5 the lattice mismatch 5 = 1.8%,® 
together with a possible misalignment of the two lattices 
given by the angle 0, create for small 9 < 2° a, periodic 
structure shown schematically in Fig. 1(a). This planar 
superlattice can be described using a set of reciprocal- 

l-(l + d)-i.R(, 




lattice vectors bn = Rnv /3 

0,1,..., 5, where R^p stands for anticlockwise rotation 
by angle ip and +0^^^ote that this set 

both rotates by (f){9) and changes its size as a function of 
9, see Fig. Bt)- which is especially important when con¬ 
sidering the zone-folding of the trigonally warped band 
structure of BLG.I^ 

In order to obtain the perturbed band structure 
and wave functions, and hence to compute the ab¬ 
sorption spectra, we use a recently proposed phe¬ 
nomenological model constructed to describe the per¬ 
turbation caused by the closely-aligned /i-BN substrate 
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on graphene electrons!^ Within this description, all 
symmetry-allowed terms containing the dominant first 
harmonics of the moire perturbation are added to the 
effective Dirac-like Hamiltonian describing electrons in 
MLG in the vicinity of one of the two corners (valleys) 
KjK' of the hexagonal Brillouin zone. This model has 
been also applied to the BLG//i-BN heterostructur^^ by 
assuming that only the bottom layer is affected by the 
/i-BN substrate. The resulting Hamiltonians 
j^BLG electrons in monol ayer and bilayer graphene, 
respectively, can be written asI 2 II 2 E 21 

iT mlg =y(j.pj.SH 

SH = (m+ f++UQf-)+ /_ + U3 /+) 

+ • [Iz X /_ -I- /+)] -I- 

f+ = J2 /- = 

n n 

pjBLG ^ I V(T p + SH + U f 

\ T1 vcr ■ p — u 

T = ^li{TzCrx-icry) + ]^V3,{ax+iT^(Jy){px+iTxPy). 

The MLG Hamiltonian ff ^^lg above is written in the ba¬ 
sis of the Bloch states on sublattices {(j){A), <p{B)} in the 
K valley and {cj)(B), —in K', whereas for the BLG 
Hamiltonian ^®lg basis is {(f){A), (j){B), (j>{A), 4>{B)} 
near K and {(j){B),—(j>(A),(f){B),—(j>{A)} near K'. In 
the BLG case, the tilde denotes the sublattice in the 
top layer. We also use two sets of Pauli matrices at, 

= (o’xjO’y), and Ti acting in the sublattice and valley 
space, respectively. The terms vcr p, v « 10®m/s,l^ arise 
due to the electron hopping from a carbon atom to one of 
the three nearest neighbors in the same layer. In MLG, 
this gives rise to the conical electronic dispersion relation. 
In BLG, due to the interlayer coupling block T, which 
includes the direct interlayer coupling 71 = 0.38eV^ 
the electronic dispersion is roughly parabolic at energies 
e <C 7 i Also present in T is the skew interlayer cou¬ 
pling, V3 « O.lrl^ which induces a trigonal distortion 
in the low energy band structure. We also include the 
interlayer asymmetry u which can be induced by apply¬ 



TABLE I: Values in meV of parameters from Eq. Q for 
the four models used in this work to calculate the absorp¬ 
tion spectra of MLG/h-BN and BLG/h-BN heterostructures. 
The given values all correspond to perfectly aligned {6 = 0) 
structures. 


model 
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Ml 

M 3 

A 

Ref. 

(1) 

11 

-21 

-18 

0 

0 

0 

0 

I12I13I 


(2) 

0 

0 

0 

-11 

-21 

-18 

0 

I12I13I 


( 3 ) 

1.26 

0.7 

-0.36 

8.98 

-7.31 

-5.63 

3.74 

117119 1 


( 4 ) 

2 

21 

-0.06 

5.2 

-42 

-5.9 

5.3 

nn 


ing electric field perpendicular to BLG and which opens 
a band gap in the electronic spectrum.l^^ISIl Finally, the 
moire perturbation SH consists of four terms: the first 
describes a simple potential modulation, the second the 
local A-B sublattice asymmetry due to the substrate, and 
the third the spatial modulation of hopping between the 
A and B sublattices. The final term, characterised by 
the parameter A, describes a non-zero average value of 
the sublattice asymmetry. Within each of the first three 
contributions to the moire perturbation, the first term in¬ 
side the round bracket, characterised by the dimension¬ 
less parameters uf, i = 0,1,3, describes the inversion 
symmetric part of the perturbation. Gorrespondingly, 
the second term in each round bracket, characterised by 
one of the dimensionless parameters u~, i = 0,1, 3, rep¬ 
resents the inversion asymmetric part of the perturba¬ 
tion. Note that in the case of BLG, the perturbation 
enters only in the part of the Hamiltonian describing the 
bottom layer since, because of the exponential decay of 



(b) 



FIG. 1: (a) Example of the graphene//i-BN moire superlattice 
(the lattice mismatch has been enlarged for visual purposes). 
The black dots denote the carbon atoms in graphene, whereas 
the red and green correspond to the boron and nitrogen in h- 
BN. The yellow rhombus depicts the superlattice unit cell, 
(b) Hexagonal Brillouin zone (BZ) of monolayer and bilayer 
graphene with the valley coordinate system Px,Py together 
with an enlarged view of the supercell Brillouin zone (sBZ) 
and its own symmetry points k and p . 
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the 2^2 orbital wave function with increasing distanc^^, 
we assume that the effect of h-BN on the top layer of 
graphene, twice as far from the substrate as the bottom 
one, can be neglected.l^ 

The position-dependent moire perturbation terms we 
wrote in Eq. Q represent all symmetry-allowed terms 
generated by the first harmonic of the moire within the 
continuum model. Various physical effects, like strain 
generated due to mutual relaxation of graphene and /i-BN 
lattice^^or electron-electron interactiorP^, contribute to 
the currently unknown values of the parameters {uf, 
u~} and A. Here, we use values as provided by four 
models published previously. All the sets, summarised 
in Table |lj were originally obtained for a MLG/Zi-BN 
heterostructu re wi th no angular misalignment {0 = 0 ). 
The first two, h^h3 | and (2), result from representing 
h-BN as a lattice of point-charges. Model (1) assumes 
that the inversion-symmetric part of the perturbation 
is dominant, whereas for model ( 2 ) only the inversion- 
asymmetric part is non-zero. The strength of the per¬ 
turbation is then governed by a single parameter, which 
we set by requiring that = 50me V for 

both of the models. Parameters for the model 
have been obtained using DFT. Model (4borrows the 
information on the stacking-dependent adhesion energies 
from DFT and uses them as input in an analytic descrip¬ 
tion of the low-energy band structure taking into account 
the elastic energy of graphene.l^ We point out that both 
models (3) and (4) explicitly include the influence on the 
electrons of deformations of the graphene lattice arising 
due to its structural relaxation on top of h-BN. Such re¬ 
laxation minimises the total energy of the crystal through 
the interplay between lowering local van der Waals in¬ 
teraction energy but increasing the elastic energy by in¬ 
troducing local strains. The evidence of such strains in 
graphene was found in Ref.|2Sl Note that it is known that 
the dominant effect of deformations in graphene crys¬ 
tals is the generation of a gauge fielcPH with sublattice 
and valley structure identical to the third term in SH in 
Eq. Q. Indeed, according to Ref. and model (4), the 
inhomogeneous, moire-periodic strain in graphene lattice 
generates significant parameters uf and u^, see Table 
|Tj This is however not the case according to model (3). 
Deformations also induce terms in SH with higher har¬ 
monics of the moir. However, these only lead to second 
order corrections to the energy of the lowest minibands. 

In order to plot the optical absorption spectra shown 
in this work, we calculate the absorption coefficient g{u}) 
for incident light with frequency uj and polarisation e = 


Snh 
g{uj) = — 


fpX' fpX 


cujS ^ hw 

P,a,A' 


+ CpA - epA' + ig 




r] is the broadening of the energy states (which, unless 
stated otherwise, we take to be ?? = 3meV). The indices A 
and A' represent the remaining quantum numbers: spin, 
valley, and miniband index. We also numerically com¬ 
pute the eigenstates of the Hamiltonians in Eq. 0 , in 
order to calculate the matrix elements of the current op¬ 
erator 

Map = (pA|j*|pA')(pA'|j>|pA), 

where ja = dH/dpa = evua are the Dirac current oper¬ 
ators. In the absence of a moire perturbation, the mono- 
layer absorption is flat with a value of Tre^ /hc W^ wh ile 
the bilayer shows a pronounced peak at huj = 71 .HHIII] 

In Fig.[^ we show optical absorption spectra calculated 
for all the four perturbation models in Table I for both 
MLG (top row) and BLG (bottom row). For the spectra 
in the left column, we assume that the structure is neutral 
(chemical potential p is at the Dirac point for MLG and 
the neutrality point for BLG). The spectra in the right 
column are for p-doped structures with hole density Uq = 
equivalent to four holes per moire unit cell (uq « 
2.5 X 10^^cm“^ for 9 — 0). If the first valence miniband 
is isolated from the rest of the spectrum (for example, 
by a secondary Dirac point somew here at the edge of the 
sBZ as suggested by experiments, I^ESl) this hole density 
corresponds to the full filling of this miniband. As a 
result, optical transitions at low frequencies are sensitive 
to the details of the reconstructed miniband spectrum. 

For the neutral structures, one expects the absorption 
spectra to be modified in the frequency range huJuLG ~ 
hvb for MLG and hiOBhC ~ for BLG, as these en¬ 

ergies correspond to optical transitions from the recon¬ 
structed electronic states at the boundary of sBZ in the 
valence band to matching states on the conduction-band 
side. However, the first three models show only weak 
modulation of the absorption spectra, both for MLG and 



where a,/3 = x,y, Cpx stands for the miniband energy 
found by diagonalisation of the appropriate Hamiltonian 
in Eq. 0 . fpX are the occupation numbers, S is the nor¬ 
malisation area of the miniband plane wave states and 


FIG. 2: Absorption spectra for MLG and BLG on h-BN with 

,2 pi 

hole density n = 0 or n = no = equivalent to filling of 

one separated miniband with holes. For clarity, the consecu- 

2 

tive spectra have been shifted by 3^^ along the y axis. 
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BLG. As a result, despite significant differences in the 
weighting of the parameters and u~ in each of the 
sets and important differences in the underlying physics 
(e.g. presence or lack of inversion symmetry), the pre¬ 
sented spectral features do not offer a clear method of 
distinguishing between the perturbation models. For ex¬ 
ample, in a recent experimental work,l^ the optical con¬ 
ductivity data was fitted using model (1) and the results 
interpreted as a confirmation of the importance of the 
sublattice mixing term, in particular u^. However, our 
results indicate that distinction between the influences 
of uf and or, as an alternative example, between 
models (1) and (2), requires more careful analysis. In 
fact, potentially very different combinations of the moire 
perturbation parameters might lead to similar spectra 
for particular choices of chemical potential. However, if 
model (4) is the physically relevant description, then the 
presence of strong absorption peaks at fiw « O.SeV for the 
monolayer, and Hu « O.leV, 0.5eV for the bilayer will be 
experimentally distinguishable. This is also evidenced 
in the right column in Fig. where we model spectra 
of MLG and BLG structures p-doped with hole density 
no- Again, the optical spectra for models (1) and (2) 
look very similar, with the most significant feature being 
a low-energy peak indicating optical transitions between 
the second and first miniband in the valence band. This 
peak, although higher and shifted to even lower frequen¬ 
cies, appears also in the spectrum for model (3). Again, 
model (4) yields the most distinctive spectrum. 

We now move beyond a “single spectrum” picture, 
which concentrated on analysing the absorption spec¬ 
trum for a given carrier density, and discuss changes in 
optical spectra as a function of increasing hole doping 
n « rig. In Fig. we show the evolution of optical ab¬ 
sorption for the four moire perturbation models as the 
doping is varied from O.Sng to 1.2no. For this density 
range, absorption spectra at low energies are dominated 
by features due to the moire perturbation (for an un¬ 
perturbed MLG/BLG spectrum with chemical potential 
/i, optical transitions with energy e < 2/j are forbidden 
due to Pauli blocking). For MLG for both models (1) 
and (2) [panels (a) and (b)], the low-energy peak due 
to transitions connecting states at the boundary of the 
sBZ in the first and second valence miniband grows in 
size with n increasing up to n ^ l.lng. A shoulder 
also develops for n Ki tiq. However, the spectrum for 
model (1) displays a clear side feature to the right of the 
main peak which appears when filling the second mini¬ 
band which is a clear experimental indicator of strong 
inversion-symmetric perturbation. For model (3), the 
peak decreases rapidly after the hole density increases 
beyond ng and splits into a double-peak feature. In 
comparison, for model (4), the height of the low-energy 
peak does not change much as a function of increasing 
hole density. Its width increases as the first miniband is 
fully filled and then decreases to roughly the same shape. 
However, similarly to model (I), a second smaller peak 
appears at slightly higher frequencies. 


For the BLG spectra, we notice a peak at lower fre¬ 
quencies than for monolayer, at ~ 50meV. For model 
(1), this peak grows for hole densities lower than rig and 
then decreases for n > uq. This is in contrast to the 
behaviour for model (2), for which this peak remains of 
similar height for n = O.Sng and n = 0.9ng but increases 
in the vicinity of ng and remains high with further in¬ 
crease of n. This is a second clear way to experimentally 
distinguish model (I) from model (2). For model (3), the 
low-energy absorption peak grows for the hole density 
increasing from n = O.Sng to n = 0.9ng but decreases al¬ 
ready for n = ng and remains small with further increase 
of the hole density. The spectrum obtained for model (4) 
is the only one for BLG that displays more clearly distin¬ 
guished features. First of all, the low-energy absorption 
peak does not change much with n < ng. However, for 
n > ng this peak splits into two, with the one at higher 
energy decaying quicker than the one at lower energy as 
the hole density is increased. Moreover, a second peak 
above O.leV exists for n < ng. 

An interesting open question is the dependence of the 
electronic moire perturbation effects on the misalignment 



FIG. 3: Absorption spectra for structures hole-doped to 

O.Sno, 0.9no, no, l.lno, and 1.2no. For clarity, the consecutive 

2 

spectra have been shifted by along the y axis. For each 
spectrum, the dashed vertical line indicates the 2/r threshold 
for Pauli blocking below which the absorption for an unper¬ 
turbed MLG/BLG would be zero. 
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FIG. 4: Absorption spectra for BLG//i-BN structures [models 

(1) and (2) only] with various angular misalignments between 

the crystalline axes of the two crystals. For clarity, the con- 

2 

secutive spectra have been shifted by 2^^ along the y axis. 


between the crystalline directions of graphene and /i-BN. 
Making use of the dependence of and u~ on 9 pre¬ 
dicted for models (1) and we show in Fig. 

the absorption spectra for the BLG//i-BN heterostruc¬ 
ture for those two models (see Ref. [5H| for the MLG/ h- 
BN heterostructure) and a range of misalignment angles 
9. Because the perturbation parameters decrease with 
increasing misalignment, the spectral features due to the 
moire weaken as well. However, the shift of the spec¬ 
tral features toward higher frequencies is visible as the 
characteristic energy of the moire, hvb, increases with 9. 
In particular, for misalignment 9 ~ 1°-1.5°, the peak at 
the energy ^ 71 due to optical transitions from the low- 
energy bands to the high-energy split band^^ is modified 
because at these angles ^ « 1 . 

One of the unique features of BLG is the possibility of 
opening a bandgap in the electronic spectrum by break¬ 
ing layer symmetry with external elec tric f ield applied 
perpendicular to the graphene layers.^^^^^ Band gaps 
Ag » |m| » 200meV have been opened with high electric 
fields .B 2 IIII interlayer asymmetry u breaks inversion 
symmetry but preserves electron-hole symmetrjESl and 
hence, in the absence of the moire perturbation, absorp¬ 
tion in the case of asymmetry u is the same as for — u.SD 
In comparison, a generic graphene-/i-BN interaction de¬ 
scribed by the parameters in Eq. ([I ) breaks electron-hole 
symmetry and opens a gap Aq quadratic in the moire pa¬ 
rameters uf, As a result, for both the moire pertur¬ 
bation and interlayer asymmetry, the gap at the neutral¬ 
ity point is Ag « |it-|-Ao|, yielding different gaps for u and 


—u. This is seen in Fig. where we show optical absorp¬ 
tion spectra calculated for BLG/Zi-BN heterostructures 
with moire perturbation effects described by models ( 1 ) 
and (2) in Table I and in the presence of |u| = SOmeV 
and |u| = lOOmeV. At frequencies close to zero, the ab¬ 
sorption falls to zercP^, indicating presence of a gap in 
the spectrum. This gap is larger for u < 0 because for 
both models Aq < 0. Because of that difference in the 
band gap at the neutrality point, although the general 
features of the spectra for u and —u are similar, those 
for u < 0 are consistently shifted to higher frequencies. 
We suggest that, if Aq is experimentally significant in 
BLG, then comparing optical absorption spectra for op¬ 
posite interlayer asymmetries u and —u should reveal the 
sign of Aq and through the dependence of Ag on uf ^ u~ 
potentially some more information about the moire per¬ 
turbation characteristics. 

Finally, in Fig.j^we present comparison of BLG optical 
absorption spectra calculated taking into account or ne¬ 
glecting the interlayer skew coupling term parametrised 
by the velocity ug in Eq. Q. For this example, we used 
parameters from model (4) and assumed chemical poten¬ 
tial to lie between the valence and conduction bands. For 
U 3 = 0 , the electronic spectrum has circular symmetry 
and so all states with the same magnitude of momen¬ 
tum p = IpI contribute to the absorption at the same 
frequency. However, nonzero trigonal warping shifts the 
energy of an electronic state by amount depending on the 
polar angle of p. As a result, both the position and inten¬ 
sity of the spectral features change. For that reason, all 
spectra in Figs. 2|5 for the BLG//i-BN heterostructure 
were carried out with ug = O.lu. 

In summary, we have discussed the optical absorption 
spectra for MLG and BLG on closely aligned /i-BN us¬ 
ing four models developed previously to describe the ef¬ 
fect of the substrate on graphene electronic dispersion. 
We showed that, unfortunately, physically very different 
moire perturbations can lead to similar absorption spec¬ 
tra. A detailed investigation of the spectra for samples 
with (almost) completely filled first miniband will pro¬ 
vide more information about the details of the electronic 
perturbation due to moire. For BLG on /i-BN, additional 
information can be obtained by measuring absorption 
spectra in double-gated devices. However, for this sys- 
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FIG. 5: Absorption spectra for BLG//i-BN heterostructures 
[models (1) and (2) only] in the presence of the interlayer 
asymmetry u. 



FIG. 6: Absorption spectra for BLG//i-BN heterostructure 
with moire parameters according to model (4) in Table I 
without (black solid line) and including the effects of trigonal 
warping. 
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tern it is important to include the trigonal warping of the 
electronic band structure during theoretical analysis. 

During the final stages of writing of this paper, we 
became aware of the work in Ref. HSJ which calculates 
the optical conductivity for MLG on h-BN within the 
framework of model (3) in the THz regime. 
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